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Abstract. A complex notion of backward stochastic differential equation (BSDE) is pro¬ 
posed in this paper to give a probabilistic interpretation for linear first order complex partial 
differential equation (PDE). By the uniqueness and existence of regular solutions to complex 
BSDE, we deduce that there exists a unique classical solution (U(t,a:) to complex PDE and 
{U(f,x) is analytic in x for each t. Thus we extend the well known real Feynman-Kac formula 
to a complex version. It is stressed that our complex BSDE corresponds to a linear PDE without 
the second order term. 
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1 Introduction 

The Feynman-Kac formula, named after Richard Feynman and Mark Kac, establishes a link 
between PDEs and stochastic differential equations (SDEs). It offers a method of solving certain 
PDEs by simulating random paths of a stochastic process. Linear real-valued Feynman-Kac 
formula was studied early in Kac (1949, 1951) as a formula for determining the distribution of 
certain Wiener functionals. Then by the theory of BSDE, Pardoux and Peng (1992) and Peng 
(1991, 1992) generalized them to nonlinear versions. They also derived some stochastic versions 
(Pardoux and Peng 1992, Peng 1992). There are many other papers in this direction, however 
we don’t list them all here. The present paper propose a complex notion of BSDE and deduce 
a complex Feynman-Kac formula for linear first order complex PDE. 

Let us be more precise. We first introduce a one dimensional complex BSDE: 


Yl’ x = + 



y: 


t.x 


Z^,T^) dr 



T t’ x clM r , t < s < T 


r , o 


( 1 . 1 ) 


where (B s ) sg j 0 T i is a complex Brownian motion with (B s ) sg j 0 ^ its conjugate counterpart, (X t,:r ) 
is defined as 

a (r) dM r + J 7 (r)dB r , t<s<T. (1.2) 

We next want to find a triple of adapted processes {(YA, T < s < T} with values in 
C x C x C which solves uniquely (1.1). We finally show that under some analytic conditions on 
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coefficients, (U(f,x) = Y(’ x ;0 < t < T, x € C} is the unique solution of the following complex 
PDE: 


I U t (t, x ) = -g (t, x , U(t, x), cr t U x (t , x), 7 t U x (i, x)), 

ju(T,x) = /i(x), 0 < f < T. ^ j 

where U takes values in C and is analytic with respect to x for each t. Note that not like the 
real Feynman-Kac formula , BSDE (1.1) gives a probabilistic interpretation to PDEs without 
the second order term, which is essentially due to the complex Ito’s formula involving analytic 
functions. We refer to Ub0e (1987), Davis (1979) and Varopoulos (1981) for complex stochastic 
analysis. 

This paper is organized as follows. In section 2, we make some preliminaries. Section 3 
proves existence, uniqueness and regularity for the solutions of complex BSDEs. In section 4, 
we establish a link between a class of linear PDEs and complex BSDEs. 


2 Preliminaries 


Let (B f ) = (Bl,B?) t > o be a standard Brownian motion in R 2 on a probability space (D, P, P) 
and (Mt)t>o be its complex counterpart, i.e. B* = B^ + iB 2 where i = \/—l is the imaginary unit. 
Let (J~^) s >t be the augmented Brownian filtration generated by the Brownian motion (B s ) s>t 
from time t. T < 0 is a fixed time. Throughout the paper we will work within the time interval 


[0, T\. For x € R 2 , x means a column vector 


x 


counterpart in C, the set of complex numbers 

= v 


Euclid norm 


lx 1 ! 2 + 


We write x = x 1 + ix 2 as its complex 

x“ ) 

For x € R 2 or x € C, We define a common 
x 2 | 2 . For x, y € R 2 , we denote by (x, y) the scalar product of x, y. 


Let M, N be two fields of (real or complex) numbers. We denote by 'Hj r (t,T: N) the space of 


; t < s < T} s.t. E 


St I Vs 


2 ds 


< oo 


jr (t, T: N) the set of continuous and progressively measurable N-valued processes 


all ^(-progressively measurable N-valued processes {(</?, 
and by Sj 

{(V> s)\t < s <T} s.t. E [sup t < s < T |^| 2 ] < oo}. C fc (M,N),C 6 fc (M,N),C^(M,N) will denote 
respectively the set of functions of class C k from M to N , the set of those functions of class 
C k whose partial derivatives of order less than or equal to k are bounded, and the set of those 
functions of class C k which, together with all their partial derivatives of order less than or equal 
to k, grow at most like a polynomial function of the variable x at infinity. 

For a single-variable function /, f denotes its derivative, for a multi-variable function /, we 
sometimes denote by f x its partial derivative w.r.t x variable. 

We now announce a result about 2x2 matrices. 


Definition 2.1 A real 2x2 matrix A is of class Cl if and only if it has the form 

( : : b ) ^ 

where a,b £ R. 

By a direct calculation, we have 

Lemma 2.1 Class Cl is an exchangeable semigroup, i.e.: let 2 x 2 matrices A, B be of class 
Cl, then A + B, A ■ B= B ■ A, \A, (\ € RJ are all of class Cl- If A 1 exists, then A 1 also 
belongs to class Cl- 

Proof. The following calculation leads to the above results: 
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a 

b 

a 

b 

a 

b 


-b 

a 

-b 

a 

-b 

a 



c —d 
d c 




1 ( a 

a' 2 +b ' 2 l _5 


f a + c 

V b + d 

ac — bd 
ad + be 



-b-d V 

a + c )’ 
—ad — be ^ 
ac — bd , 


c —d \ { a ~b\ 
d c J \ b a J ’ 


3 Complex BSDE: existence, uniqueness and regularity 

Let X. t,x satisfy (1.2) with u(-), 7 (•) being deterministic and square-integrable functions. We 
introduce the following complex BSDE: 


r T r T 

Yi’ x = h(X%, x ) + g (r, X*’*, Y *■*, Z £ >*, T £ ’ x ) dr- Z, 

J s J s 


Yc 


r T 

J T^c 


t < s <T (3.1) 


where h : C —> C, 5 =<7 (r, x) + a (r) y + f3 (r) z + 0 (r) 7 , g: [0, T] x C — > C. 

We assume the following 

(HI). (Polynomial growth) g(s , •) € Cp(C,C) and h € Cp(C,C). 

(H2). (Bounded derivatives) a, (3,6 are C-valued bounded and deterministic functions 
(H3). (Analyticity) h,g are analytic w.r.t spatial variable, e.g. for x=x 1 + ix 2 , 

dh 1 _ _ dh 2 
dx' A dx 1 ' 

Remark 3.1 From the Maximum Modulus Principle, for a nonlinear complex function g (r, X, ¥, Z, T), 
the derivatives of g w.r.t ¥, Z,T are bounded and analytic on the whole complex plane means 
that g is linear in ¥, Z, T. 

BSDE (3.1) is equivalent to the following 2-dinrensional BSDE: for s € [0, T] 


Y? 

X 2 


or the following 


h x Xl, Xf 
h 2 (X^X 2 ) 

Zl -Z 2 


+ 


L 


T / r .l ( r yl y2 yl y 2 7 1 7 2 F 1 F 2 ! 

y y 1 y *-r 1 i 1 r "> ± r i ^r i ^r ? r ? r J 

„2 / yl X 2 V 1 V 2 7 1 7 2 F 1 r 2 l 

y y i y *-r i y *-r i 1 r i 1 r > ? 1 r? 1 r/ 


dr 


Zj 


dR 

dB 


rj r 2 
r, 2 -r 1 


dR 

dB 


(3.2) 


y 1 

s 

x 2 


, 5 , _ f d 1 {x^x 2 ) l + 


where X = 


h 2 (X r \X 2 ) 

,Y S = 


X 1 

X 2 


f 1 (r, X r , Y r , Z 
/ 2 (r, X r , Y r , Z 


dr — 


r 


y 1 

y2 


, Z = (Z 1 ,Z 2 ) = 


z 11 z 12 
z 21 z 22 


'11 

V 

z 21 N 

\ ( dBl 

'12 

V 

^ 22 , 

I ^ dB 2 

| 

z 1 

+ r 1 

r 2 + z : 

r 2 

-z 2 

z 1 -r 


(3.3) 


f 1 (r, X r , Y r , Z r 


= ^(fX^x 2 ,!; 1 ,!; 2 , ^(z r n + z 22 ), i(z r 12 - z 21 ), 

^(Z r n -Z 22 ),i(Z r 12 + Z 21 )), 


(3.4) 


/ 2 (r, X r , Y r , Z r 


= 5 2 (r,x r \x 2 ,y r \y r 2 , i(z r n + Z 22 ), i(z T 12 - z 21 ), 
i(z r u -z r 22 ),i(z r 12 + z r 21 )). 


(3.5) 
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One can check that / = (Z 1 ,/ 2 ) is analytic in (X,Y,Z), e.g. = ~^r, 

and satisfies the Lipschitz condition in (Y,Z), thus there is a unique pair ^Y^Z^^ € 
Sjr (t, T; R 2 ) xVffplt, T; R 2 xR 2 ) which solves the real BSDE (3.3) (see Pardoux and Peng 1990). 
Therefore there is a unique triple ^Y^’ x , Zs ' x , T^’(f, T; C) x T\ C) x T ; C) for 

BSDE (3.1). 

Remark 3.2 T/ie analyticity is not used for the existence and uniqueness of the solutions for 
BSDE (3.1). It is just useful when we derive the Feynman-Kac formula. 

Theorem 3.1 {Ys’ x ;t < s < T, x € C} is analytic in x and continuous in ( s,t ). 

Before proceeding to the proof, we first state a useful corollary: 

Corollary 3.1 For any t S [0, T], the mapping x —>• Y^ is analytic, the function and its partial 
derivatives of order one and two being continuous in (t,x). 

Proof of Theorem 4.1. It suffice to prove the following 
Step 1. {Y s ’ } is continuous in (s,t), { d s x } is continuous in x. 

Step 2. Cauchy-Riemann equations: = (Py, 2% = — 

Step 1 follows immediately the fact that 


E 


SU P 4 / U< s <tI Y s - Y fr < Cp(l + I*| 9 )(|x - X| p + It - t\ 2 ) 


E 


SU Pt A t< s < T l ^h Y s 


Ai Yl’ x - AiY* ,£ | p 


h i s |-j < c p (l + \x\ q + \x\ q + \h\ g + \h\ q )(\x-x\ p + \h-h\ p + \t-t\ 2 ) 

where AfYl’ 1 = (Ys’ x+he — Y t s ’ x )/h, h € R\{0}, {e 1 ,e 2 } is an orthogonal basis of R 2 . 

We now prove the Cauchy-Riemann equations. Pardoux and Peng (1992, eq. 13) says that, 

VY(- J = f T [f x (r, X£ x , Y*> x , Z(Y) VX 4 ; 1 + /' (r, X(Y, Y 4 ’*, Z*> s ) VY£> x 

J S 

+/' (r, X): r . Y'- ;r . Z‘- x ) VZ*’ x ]dr + F ((VZ 4 ’*)*, dB r ) ,se[t,T). (3.6) 

J S 

where VY^ is the matrix of first order partial derivatives of Y^’ x (x denotes the initial condition 
of SDE (2.1)). VXs’ J and VZ)Y are defined analogously. Let (M))t < s <7 be the solution of the 
following matrix-valued SDE: 


dM s = M s fyds + (M s / 4 ,dB s ) ,t < s < T. 

= I 

which is given by 



2 Jz 


c r»r 


dr + 


j* (f z ,dB r )},t<s<T. 


Since /', /' is of class Cl, Mf is also of class Cl. 


(3.7) 


(3.8) 
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Then applying a 4-dimensional version of Ito’s formula (see 0KSENDAL 2005, Th.4.2.1) to 
MjVYs’ 1 , we deduce that 

d(M s VYi’ x ) = M s d(VY‘ ,:r ) + (dM s )VY$>* + (dM s ) s d(VY 4 ’ x ) 

= —M s (/'V YY + /'VZ*-* + /'VX 4 /) ds + <M S VZ 4 /, dB s ) 

+M S /'VY 4 - X ds + (M s /'VY 4 ’ x , dB s ) + Msfyz'fda, s € [t, T\. 

Therefore 

/■T /-T 

M s VY 4 ’ i = M t VY^+ 1 M r /'vx*’ x dr- / <M r (/'VY*’ x + VZ 4 /)*, dB r > , 

Js J s 

then 


VY 4 ’ X = E 


m T s h! (x^f) VX!f x + f M^/'VX^dr | E r 

«/ S 


(3.9) 


By Lemma 2.1, It is known that VY)’ 1 is of class C l, thus the Cauchy-Riemann equations hold 
true. The proof is complete. □ 


Remark 3.3 Let Zf = Z 2 ) = ( 

Zj c Zj , 


11 712 \ (a 

21 7 22 h then f'z = ( A , B ), (/')* = 



B 

df 1 ap 
az^r az^ 
a / 2 af 2 
az 51 az 27 


Smce / is analytic in Z = Z l + i ■ Z 2 , one can check that A, B are of class Cl- 


For BSDE (3.3), Pardoux and Peng (1992) proved that: for any 0<f<s<T, xG R n , 
(Z 4 ’ ,x )* = VYs’ x (VXs X ) _1 ff* and particularly (Z^’ x )* = VY l’ x (T*. Note that VX^’ 2 = I in this 
paper, for BSDE (3.1), we have the following results. 

Proposition 3.1 For 0<t<s<T,x€C, 


dY t,x 

\^S ) — 


(Tf) = 


dx 

dY 4 / 


dx 


-7s, 


Proof. By Pardoux and Peng (1992, Lemma 2.5), we get that 


t (z 1 ^ r 1 r 2 -z 2 
J “ r 2 + z 2 z x — r 1 


From the above equation, we have that 


dY 1 

dY 1 

dx\ 

dx\ 

dY 2 

dY 2 

dx 1 

dx ' 2 


a 1 +yi y 2 — a 2 
7 2 + a 2 a 1 — 7 1 


(3.10) 

(3.11) 


(3.12) 


Zl + ^ “ ^ (CJ1 + 7 ' } + ^ {a2+ 72) ’ 


dY 


dY 1 


z ‘- r ‘ = ^'-D + ^ 2 -D. 


(3.13) 
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and 


Therefore, 


and 



r 2 + z 2 


dY 2 
dx 1 
dY 2 
dx 1 


(7 1 - cr 1 ) + 

(7 1 + CT 1 ) + 


dY 2 
dx 2 
dY 2 
dx 2 


(7 2 -A 

(7 2 + a 2 ). 


Z l 

r 1 


dY 1 , aY 1 2 
ax 1 ^ dx 2 a ’ 
aY 1 j aY 1 2 

ay 7 + ay 7 ’ 


(3.14) 


(3.15) 


Thus 



dY 2 , aY 2 2 
ax 1 ^ 9x 2 17 ’ 

ay 2 j ay 2 2 
ay 7 + ay 7 ' 


(3.16) 


Z 1 + * ■ Z 2 = ( 


dY 1 

dx 1 


o - 1 + 


dY 1 

dx 2 


(j 2 ) + i{ 


dY 2 

dx 1 


a 1 + 


dY 2 

dx 2 


.dY 1 .dY 2 .. , . 2 \ 

<r ‘ J = ( a?- + ! a ?' )(<7 +!<T) ' 


R+i 


v a ® 1 


that is 


, dY 1 .ay 2 , . 2 \ 

( ay + *ay )(7 + * 7 } ’ 


(zl**) 

(Tf) 


dY *’ 1 

dx 

dY*’* 


□ 

By Proposition 3.1 we know that, if <r = 0(resp. 7 = 0 ), then Z = 0(resp. T = 0). Since 
there is a unique solution (Y, Z,T) for BSDE (3.1), we have the following results: 


Corollary 3.2 Under conditions (HI) ~(H3), there is a unique solution in Sy (f, T;R 2 ) x 
T-Ljrit, T;R 2 x R 2 ) respectively for the following two real forward-backward SDEs: 


X] \ 


1 2 
<Jr —^r 

\ ( dB) \ 


Xl ) 

2 1 

Cj. a}; 

) V dB 2 r ) ’ 


Yl ' 

\ ( id {x),xl) 

H( 

„1 ( r yl y2 yl y2 y\ 
y y ? i 1 r > J r 5 ^r 5 

y2 pi p2\ \ 

ZJ T , X r , 1 r j \ 

zl,rj.,r 2 r ) ) 

. V? , 

) V h 2 (x),xl) 

n 2 ( r yl y2 yl y2 yl 

y ^5 A r , ij, , ij. , ZJ r , 


-Is 


Z) -Z( 

zl zl 


dBl 
dBl 


s € [t, T] , 


s € [t, T] , 


(3.17) 
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and 


X] 

X 2 s 

y} 

n 2 


, r s ( T'r Tr 
+ jt v 7 r 
h 1 {X},X?) \ , 

h 2 (x),X?) ) + X 

r T ( rl r 

Js U? 


r 2 

i r 

-r 1 

1 r- 


dBj. 

dB 2 

n l / r v-l y 2 V 1 V 2 7 1 

y y ') yv r'>- /v r'> 1 r i 1 r > ^r 
n 2 ( r yl v2 yl y2 yl y2 
y \' j - /v r 5 -y*-r i 1 r ? 1 r i ^r ? ^r 


^2 t - ' 1 p2\ 
r > ’ 1 r’ i r; 

r 1 r 2 ) 

X r , X r I 


dB). 

dBl 


dr 


s £ [t, T ] , 


5 € [7, T] , 


(3.18) 


Remark 3.4 The above two BSDEs are real-valued BSDEs with Z-constraints. Usually there 
are no solutions for constrained BSDEs. 


4 Complex BSDE and associated PDE 

Consider the following complex PDE: 

f U t (f, x) = -g (t, x, U(t, x),a t V x (t, x), 7 ' t U x (t, x )), 
U(T, x) = h(x), 0 < t < T. 

where U : R + x C — > C, a = a 1 + ia 2 , 7 = 7 1 + ry 2 , cr 1 7 1 = a 2 7 s . 


(4.1) 


Theorem 4.1 Let h,g satisfy (HI) ~(H3). If for each (t,x), {U (t,x);t < s < T,x € C} is 
analytic with respect to x and continuous in t, and satisfies PDE (f.l), then 

(Yi-xfXn ■■= {v(sX’ x )^sV x (sxxusVxX’ x )) 

solves BSDE (3.1). Furthermore, 

(v(t,x),a t is x (t,x),”/tVx(t,x)) 0 < t < T = (yYXfXf 

and 


U(t,x) =Y\’ X = E 


h (X^) + J g (r, X£*, Y(.’ x , Z*’®, T£>*) dr 


Before proving the above theorem, we need an Ito’s lemma. 

Lemma 4.1 Let X* € S(p (0, T;C), bt, atilt x Hjr(0, T; C), such that cr 1 7 1 = a 2 7 s for almost 
all t and 

dKt = btdt + atdM t + 7 tdM t , t > 0. (4-2) 

If¥(t,x) = u(t, x 1 , x 2 ) + iv(t, x 1 , x 2 ) is an analytic function w.r.t the complex variable x and 
continuous in t, i.e. for any t, Vx € €, F satisfies the Cauchy-Riemann equations: 

du _ dv 

dx 2 


= then 


d¥ <9F 

dF(t,X t ) = —(t,Xt)dt + —(t,X t )dKt 


d¥ <9F <9F d¥ 

= —(t,X t )dt + —(t,X t )btdt + —(t,Xt)a t dMt + -^(t,X t )i t dM t ,t>0. (4.3) 

where ||(f, x) is the complex partial derivative of F w.r.t x. 





Proof. 


dF(t,X t ) = du(t,Xl,X?)+idv(t,X},X?) 


du du o 9ii , 1 . 

a?''* 1 + + W* + 2 ( a7r 


<9 2 m , <9 2 u 


d < .X'l >i + 


+* 


<9u <9u 9 du , 1, <9 2 u 

a^ dx ‘ + + M dt + 2 ( a(x I p 


L f 

<9u 


d < X 1 > t 4 


0(x 2 ) 2 

d 2 v 


d<X 2 > t ) 


d(x 2 ) 2 


d < X 2 > t ) 


/ du 

ycte 1 ^cte 1 


^ (dXt + idX 2 ) + ^ 


/ <9u dv\ 1 

m + ‘m' dt 


<9F <9F 

— (t,X t )dt + —(t,X t )dX t 


where we have used the conjugate harmonicity of function u, v and the condition cr 1 7 1 = 
er 2 7 2 . □ 

Proof of Theorem 4.1. It suffices to show that 


{iJ(sX’ X )^sV x (sX’ X )^sV x (s,X t s ’ x y,t< s < T} 

solves BSDE (3.1). Applying the complex Ito formula to U(s,Xs’ a: ) between s = t and s = T, 
we get that 


¥*/ = h (x^) + f T [g (r, X^, U(r, X^), a r V x (r, X^), 7r U x (r, X^))] dr 


<r r U x (r, X: 


t,x\ 

y J UjlLDj' 


7 r .U a; (r, X(. ,3: )(i]B r , t < s <T. 


Thus (¥*’ a ,Z FX*) = {U(s,X t /),cr s U a; (s,X t /), 7 s U x (s,X t s ’ a; )} solves BSDE (3.1). □ 
We now show the converse of Theorem 4.1. 


Theorem 4.2 Let h,g satisfy (HI) ~(H3). Let {(¥s ,a: );f < s < T } be the solution of BSDE 
(3.1). Then (U(t, x)) 0<t<T = is the unique classical solution of backward PDE 

(4-1) and U(t, x) is analytic in x for each t. 

Proof. Uniqueness follows from Theorem 4.1. We now prove that is a solution to 

PDE (4.1). Let 5 > 0 s.t. t + 5 < T. Clearly Y^ = ¥ ( ^’ t+s . Hence by the complex Ito 
formula and the analyticity of U(t,x) in x and BSDE (3.1), we have 


U(t + 5, x) — U(t, x) = 


\J(t + 5,x)-V{t + 6,Xl'+s) + U(f+ <5,X^) - U(t,s) 

rt-\-5 

- g(r,x, X*’ 2 , Y*’ 2 , Z f r ’ x , T*’ 2 ) dr 


rt+S 


a r V x (t + 6, X; 


rt+S 


7rU x (t + s,: 


rt.x 


rt+S rt+S 

+ / Zl’ x dE r + / T ? r ’ x c 

Jt Jt 
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Let t = to < t\ < ... < t n = T, we get 

ra-l r t 


h(x)-U(t,x) = / g(r,xXr X Xr X ,K ,X X- X )dr 

i =0 ' t% 

ru+i 


J- 

E/ 

j=0 

+ [ U+1 [Z^ - a r V x {ti+i,K’ x )] 

Ju 

+ [ 1+1 [T^- 7r U^ m ,X^)] 

Ju 


Let the mesh size sup 0 <j< n _i (L+i — tf) —>• 0, we obtain the limit 


U (t,x) = h(x)+ / g (r, x,U(r, x), cr r U x (r, x), r y r ^x(i r ', %)) dr, t<s<T. 


(4.4) 


□ 


Remark 4.1 Viscosity solution for PDE (4.1) is not involved in the present paper because, 
the analyticity of parameters of PDE (4-1) leads to existence of the first and the second order 
derivatives. 


The well-known Cauchy-Kovalevski theorem states a local existence and uniqueness of solu¬ 
tion for partial differential equations whose coefficients are analytic functions, associated with 
Cauchy initial value problems. A special case was proven by Cauchy in 1842, and the full result 
by Kowalevski (1875). Theorem 4.2 extends the first order Cauchy-Kovalevski theorem to the 
case of global solutions. 
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